Draft version March 2, 2013 

Preprint typeset using I^T^X style cmulatcapj v. 5/2/11 



DEEP ABSORPTION LINE STUDIES OF QUIESCENT GALAXIES AT Z ~ 2: THE DYNAMICAL MASS-SIZE 
RELATION, AND FIRST CONSTRAINTS ON THE FUNDAMENTAL PLANE 

S. Toft 1 , A. Gallazzi 1 , A. Zirm 1 , M. Wold 1 , S. Zibetti 1 - 2 , C. Grillo 1 , A. Man 1 

Draft version March 2, 2013 

ABSTRACT 

We present dynamical and structural scaling relations of quiescent galaxies at z = 2, including the 
dynamical mass-size relation and the first constraints on the fundamental plane (FP). The backbone 
of the analysis is a new, very deep VLT/X-shooter spectrum of a massive, compact, quiescent galaxy 
at z = 2.0389. We detect the continuum between 3700-22000A and several strong absorption features 
(Balmer series, Ca H+K, G-band), from which we derive a stellar velocity dispersion of 318 ± 53 km 
s . We perform detailed modeling of the continuum emission and line indices and derive strong 
simultaneous constraints on the age, metallicity, and stellar mass. The galaxy is a dusty (Ay — 
0.771q 32) solar metallicity (log(Z/Z Q ) = 0.02^q'4°) post starburst galaxy, with a mean luminosity 
weighted log(age/yr) of 8.9 ± 0.1. The galaxy formed the majority of its stars at z > 3 and currently 
has little or no ongoing star formation. We compile a sample of three other z ~ 2 quiescent galaxies 
with measured velocity dispersions, two of which are also post starburst like. Their dynamical mass- 
size relation is offset significantly less than the stellar mass-size relation from the local early type 
relations, which we attribute to a lower central dark matter fraction. Recent cosmological merger 
simulations qualitatively agree with the data, but can not fully account for the evolution in the dark 
matter fraction. The z ~ 2 FP requires additional evolution beyond passive stellar aging, to be in 
agreement with the local FP. The structural evolution predicted by the cosmological simulations is 
insufficient, suggesting that additional, possibly non-homologous structural evolution is needed. 



1. INTRODUCTION 



Since first being disco vered in NIR surveys (e.g., 
FIRES Franx et al. 2003), massive, evolved galaxies at 
z ~ 2 have been a subject of extensive research, char- 
acterizing their properties fr om increasingly larger and 



more complete samples (e.g. |van Dokkum et al. 2006 
Franx et al.|[2008{ |ToTt et al.| \IW>\ |McCracken et al 



20101 [Williams et al.||2010| |Newman et al.||20111|Dam- 
janov et al.J2011||Cassata et al |2011[|Wuyts l et al.|201ip 



By combining photometric redshifts, stellar masses, and 
star formation histories derived from Spectral Energy 
Distribution (SED) fits to broad band photometry in- 
vestigators have probed their general properties: They 
are massive (M> 10 11 M fl ) galaxies, with space densities 



5 x 10~ 4 Mp c~ 3 (|vaii Dokkum et al_ 
et al.||2009|) which cluster strongly (Daddi et al. 



2010 Marchesini 



Quadri et al.|20U8|) and approximately 30 — 50% o: 
are q uiescent, with little or no star formation (e.g 



2003 
them 



et al. 2006 Toft et al. 2007 2009 Kriek et al. 



Kriek 



2009). 



In high resolution imaging, the star forming galaxies 
have extended and in some cases disturbed morpholo- 
gies, while the quiescent galaxies are uniformly extremely 
compact (effective radius r e < lkpc), corresponding to 
average stellar mass densities (within r e ) 1-2 orders of 



magnitudes higher than in local early types flDaddi et al 



20051 IToft et al 



20071 (Toft et~aTT 



2005; ITrujillo et al 



2007T|Cimatti et al.| 



2006 



2008 



20091 |Conselice et al.|)2011[ |Cassata et al.||2011[ [Wuyts 



Zirm et al. 
rick et al. 
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et al.|201lb. but see also |Mancini et al.j |2010| ) ; |Newman 
et al.[ ( |2010[ ). Furthermore, there is evidence that a sig- 



nificant traction of z ~ 2 quiescent galaxies have fiat, 
exponential disk like surface brightness profiles, rather 
than the r 1 / 4 nrnfiles found in local nniescent ca.laxies 



(Toft et al. 2005, 2007; 


|van Dokkum et al.||2008 


Wei et al. 2TJTT 


"| Wuyts 


et al.||201l|) 



van der 



If galaxies with these proper ties exist in the local uni 



verse, they are extremely rare jTrjijillo et all|2009l |TaF 
lor et al.|20i0] |Shih fe Stockton|2011 |). The evolutionary 



path ot the compact quiescent galaxies to the local uni 
verse is not well understood. Their high stellar masses, 
compact morphologies and quiescent nature makes it nat- 
ural to assume they are the progenitors, or "seeds" of 
massive local early type galaxie s, which is also supported 
by number density arguments ( van Dokkum et aL|2010 1 
but they have to go through significant structural evo- 
lution, increasing their effective radii by factors of 3-6, 
and possibly transforming their flat exponential disk like 
structure into bulge-dominated/early type systems. In 
the following we will refer to these massive quiescent 
compact z ~ 2 galaxies as SEEDs. Minor and major dry 
merging are the primary candidate processes for puffing 
up their sizes and transforming their profiles from flat 



disk-like systems to spheroidal early types (e.g. Khoch 



far fe Silk 2006] |Naab et al |2009| [Bezanson et al.|2009 ) 



One of the mam uncertainties m the interpretation of 
the properties of SEEDs is that it is mainly based on pho- 
tometric redshifts and SED fits to broadband photome- 
try Due to their quiescence, they are void of strong emis- 
sion lines and have very faint rest-frame ultraviolet (UV) 
continua. It is challenging to detect absorption lines, that 
can be used for redshift determination of quiescent sys- 
tems. The red rest-frame optical color of SEEDs makes 



2 
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them brighter in the observed near-infrared (NIR), but 
due to the limited sensitivity of NIR spectrographs, spec- 
troscopic confirmation have only been possible for small 
samples of the very brightest examples. In a NIR spectro- 
scopic survey of nine of the brightest examples, using the 
Gemini Near Infrared Spectrograph (GNIRS), |Kriek et| 
al. ( 2006 1 were able to detect the continuum well enough 
to determine the redshift from the position of the 4000 A 
break, and confirm the relatively old age (1-2 Gyr) and 
quiescence from the shape of the continuum and lack of 
strong emission lines. Even with relatively long exposure 
times (1-4 hours) on an 8m telescope it was not possible 
to make significant detections of absorption lines. Only 
in a much deeper (29 hour) follow-up spectrum were ab - 
sorption lines detected for one object ( Kriek et al.|2 009). 
In this case, it was also possible to measure the stellar 
velocity dispersion of a z ~ 2 galaxy for the first time. 
The dispersion is very high, consistent with what was 
expected from its high derived stellar mass density ( Toft 



et al.|2007||Cimatti et al.|2008|[van Dokkum et al.|2009p 

Since then velocity dispersions have been derived for a 
few handful of quiescent galaxies in the redshift range 
1 < z < 2 dCenarro & Trujillo||2009| |Cappellari et al. 



20091 INewman et al.p010| |Onodera et al.p010| |van de 
Sande et al.||201ip . 



With the advent of new, more sensitive spectrographs 
it is now feasible to obtain "absorption line quality" spec- 
tra for samples of SEED galaxies in reasonable amounts 
of time. In this paper we present the first results of such 
a survey that we conducted with the X-shooter spectro- 
graph on the Very Large Telescope. As this is the first 
detailed spectroscopic investigation of this type of galaxy, 
we describe stellar population and velocity dispersion fits 
in considerable detail. 

The structure of the paper is as follows. In Section [2] 
we describe the target selection, acquisition and reduc- 
tion of a deep X-shooter spectrum of a quiescent galaxy 
at z ~ 2. In Section [3] and [4] we perform a detailed 
analysis of the spectrum, fitting the UV-NIR continuum 
and absorption lines indices, with stellar population syn- 
thesis models, to derive stringent constraints on the age, 
stellar mass, dust content and metallicity of the galaxy. 
In Section [S] we fit the surface brightness profile of the 
galaxy to derive its structural parameters. In Section [6] 
we measure the velocity dispersion from the absorption 
lines, and derive the dynamical mass. In Section [7] we 
compile a sample of 3 other z ~ 2 galaxies from the lit- 
erature with measured velocity dispersions, and compare 
scaling relations between their structural and dynamical 
properties with those observed in local galaxies, to derive 
constraints on the dynamical-mass size relation and, for 
the first time, the fundamental plane relations at z=2. 
Finally we summarize and discuss the results in Section 
[5J Throughout the paper we assume a standard flat cos- 
mology (0 A ,fi M = (0.7, 0.3), # = 73 km s^Mpc 1 



All m agnitudes are referenced to the AB system (Oke 
19741 



2. OBSERVATIONS 

2.1. Target selection 

Our objective was to take advantage of the unique com- 
bination of sensitivity, spectral coverage and resolution 




K-band image 



r neighbor model 



Fig. 1.— Left: K-band image ( 10.7" x 10.7") of UDS19627 (WF- 
CAM), with the position of the X-shooter 11" X 0.9" slit indicated. 
North is up, East is to the left. The seeing is 0.7". The brighter 
galaxy north-west of the target is at lower redshift. Middle: We 
modelled the surface brightness profiles of the two galaxies simul- 
taneously with GALFIT. In this image the best fitting model of the 
neighbor has been subtracted, demonstrating its minimal influence 
on our fit of the main target. Right: Residual after subtracting 
both best fitting surface brightness models. The red cross indi- 
cates the centroid of a 24^m source associated with the foreground 
galaxy. 

of X-shooter to obtain a spectrum of a SEED galaxy 
of unprecedented quality, to allow a detailed study of 
the absorption lines for the first time, and to obtain a 
spectrum of sufficient S/N to be used as a template for 
designing programs to observe larger samples. 

The target UDS19627 (RA=02:18:17.1, DEC=- 
05:21:38.8 ) was chosen from the UKIRT Ultra Deep Sur- 
vey (UDS), as a (K-band) bright examp le of a compact, 
quies cent, massive galaxy with z ~ 2 (Williams et al. 
20101. Based on the broad band photometry (See Ta- 
ble :?), the target is a very massive, quiescent galaxy 



at z 



phot 



= 2.021^ with K=20.19±0.01. Images of the 
target are shown in Figure [l] The bright galaxy ~ 3" 
north-west of the target is at lower redshift (UDS19771, 
Zphot = 0.58, Wil liams et al.|[2010 ). This galaxy may 
be boosting the brightness ot our target slig htly , through 



the gravitational lensing effect (see Section 7. 



2.2. X-shooter Data 

We obtained deep spectra with the X-shooter spec- 
trograph, mounted on ESO's Very Large Telescope. X- 
shooter is a single object echelle spectrograph, consist- 
ing of three arms which cover simultaneously the spec- 



tral range 3000-25000 A ( |Vernet et al~]|201l[ ). We ob- 
tained a total exposure time of five hours over three 
nights from October 19-21 2009, as part of the X-Shooter 
GTO program. Details of the observations are given in 
Table [TJ On the first night the seeing was very good 
(FWHM~ 0.5"), so we used a narrow slit. On the fol- 
lowing nights the seeing was 0.7" — 1.0", so we switched to 
a larger slit width, to minimize slit loss. We performed 
a blind offset to the object from a nearby bright star, 
and nodded on the slit with a ABA'B' sequence with 4 
independent positions. 

2.3. Reduction 

We reduced the data using a combination of the ESO 
X-shooter pipeline (version 1.2.0 in physical mode) and 
our own scripts. We experimented with different reduc- 
tion recipes. The cleanest, highest S/N spectrum (be- 
tween the skylines) was achieved when reducing the in- 
dividual exposures in STARE mode. To determine accu- 
rate offset between the trace in different exposures, we 
"collapsed" regions between skylines, in the rectified, sky 
subtracted frames and measured the shifts between the 
peaks of the spectral point spread functions (SPSF, the 
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TABLE 1 

Details of the X-shooter observations 



Night 


Exptimc [s] 


Slit size["] 






[UVB/OPT/NIR] 


[UVB/OPT/NIR] 




10/19/2010 


4 x 800/4 x 800/4 x 900 


0.8 x 11/0.9 x 11/0.6 x 


11 


10/20/2010 


8 x 800/8 x 800/8 x 900 


1.0 x 11/0.9 x 11/0.9 x 


11 


10/21/2010 


8 x 800/8 x 800/8 x 900 


1.0 x 11/0.9 x 11/0.9 x 


11 



TABLE 2 

Photome try: A B magnitudes for UDS19627 from the UDS survey QWilliams et al.| 
|2010[ ) . ChI and Ch2 refer to Spitzer/IRAC 3.6 and 4.5/xm bands 



u 


B 


V 


R 


I 


z 


J 


H 


K 


Chi 


Ch2 


25.72 
±0.21 


24.76 
±0.03 


23.94 
±0.03 


23.47 
±0.03 


23.15 
±0.03 


22.62 
±0.03 


20.94 
±0.03 


20.46 
±0.03 


20.19 
±0.02 


19.67 
±0.03 


19.58 
±0.04 



object is effectively a point source due to the seeing). We 
then shifted the individual exposures to the same posi- 
tion and combined them by computing the sigma clipped 
average of each pixel (masking out bad pixels and cosmic 
rays). We also created a 2D variance map in which each 
pixel is the variance of the 20 pixels in the stack. 

The individual background subtracted, rectified 2D 
NIR frames were corrected for telluric absorption, using 
telluric star observations that were observed immediately 
before or after each 1 hour observation block. 

This correction was then applied to each row in the in- 
dividual 2D frames before they were combined. Finally, 
we extracted the ID spectrum, using optimal extraction 
(Home 1986), weighing the pixels both by the SPSF 
and by their variance. For flux calibration we used a 
spectrum of the spectro-photometric stand ard standard 
BD +17°C4708 ( |Bohlin fe Gilliland||2004) in combina - 
tion with the CALSPFC HST database ( |Bohlin||2007| ). 
To correct for slit losses we compared to the broad band 
photometry, and applied a constant scaling to each of the 
three arms. 



3. SPECTRAL ANALYSIS 

In Figure [| we show the full composite 3500-20000 A 
X-shooter spectrum, constructed by stitching together 
the flux calibrated UVB, VIS and NIR arm spectra. The 
spectrum confirms the expectation from the broad band 
photometry, that the galaxy is a massive quiescent galaxy 
at z ~ 2, with a strong 4000 A break in the J-band. 
Interestingly the (flat) restframe UV continuum of the 
galaxy is detected throughout the VIS and UVB arm, 
all the way down to ~ 3700 A. The overall shape of the 
continuum is characteristic of a post starburst or E+A 



galaxy (e.g. Quintero et aLp004[ ) rather than of an early 
type galaxy! A number of strong absorption lines are 
detected (Balmer series, Ca H+K, G-band). From the 
absorption lines we derive a redshift of 2.0389±0.0004. 

4. STELLAR POPULATIONS PROPERTIES 

In order to constrain the galaxy stellar populations, 
we analyse the X-shooter spectrum by means of a Monte 
Carlo library of 1 00000 model spectra, based on |BruzuaI| 
fe Charlot| ( [2003l hereafter BC03) Simple Stellar Popula- 
tion (SSF) models co nvolved with ran dom star formation 
histories. We adopt a |Chabrier ( |2003 ) Initial Mass Func- 



tion (IMF). The SFHs are modeled with an exponentially 
declining function on top of whic h random bursts of sta r 
formation can occur. Following |Charlot fc Fail] ( |2000|) 



attenuation by dust is modeled by means of the total ef- 
fective optical depth Ty experienced by young (< 10 7 yr) 
stars still in their birth cloud and the fraction [i of it 
contributed by the ISVQ 

We adopt a Bayesian statistical approach in which 
all models in the library are compared to the observed 
galaxy. The likelihood of each model is computed as 
exp(— x 2 /2) by comparing a selected set of observables 
with the model predictions. The posterior probability 
density function (PDF) of each physical parameter of in- 
terest is computed by weighing all the models by their 



ters (see e.g. 


Kauffmann et al.|2003 




Gallazzi et al.|2005 


Salim et al. 


2005). The estimatec 


parameter is taken 



confidence interval is estimated from the 16 th and 84 th 
percentiles of the PDF. In particular we derived PDFs 
of the luminosity- and mass-weighted ages, the stellar 
metallicity, the stellar mass and the parameters describ- 
ing dust attenuation. 

In order to compare the data with the models we 
use a lower resolution version of the X-shooter spec- 
trum binned over 18 pixels which corresponds to a 
rest-frame AA=3A comparable to the spectral resolution 
of the BC03 models in the optical range. We adopt a 



3 We adopt the model library used is |Salim et al.| | |2005| , but 
we restrict the formation time to vary between the age of the Uni- 
verse at the redshift of the galaxy (i.e. 3.2 Gyr) and 1.5 Gyr in 
order to reproduce the current quiescent nature of the galaxy as 
witnessed by the lack of emission lines (we checked that the precise 
choice of minimum formation time does not affect the stellar mass 
estimates and only minimally the age estimates). The other param- 
eters defining the SFH and dust att enuation follow the same prior 
distributions as in |Salim et al.||2005| ). The star formation timescale 
can vary between and 1 Gyr -1 . The probability of having a burst 
is set such that 50% of the models experience a burst in the last 
2 Gyr. Bursts are parametrized in terms of the fraction of stellar 
mass produced, which is logarithmically distributed between 0.03 
and 4, and their duration can vary between 3 X 10 7 and 3 X 10 s yr. 
The metallicity of each model can vary between 10% and 2 times 
solar and it is fixed along the SFH. The total effective optical depth 
Ty can vary between and 6, and /i can vary between 0.1 and 1. 

4 The binned spectrum is computed as the noise-weighted me- 
dian flux in each bin and the error is obtained from the noise- 
weighted 16 th and 84" 1 percentiles. 
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Fig. 2.— Full 2D (top) and extracted ID (bottom) X-shooter UVB (blue), VIS (green) and NIR (red) spectrum of UDS19627 (corrected 
for slit losses). The colored curves show the full resolution (0.5 A) spectrum, the black curve is binned to a resolution of 10 A, where the 
value in each bin is the weighted median of the indiv idu al pixel values. Also plotted are the broad band fluxes (blue circles). A number of 
absorption lines are detected (see also Figures [9| and [lOp 



larger AA of 16A in the rest-frame UV (UVB and VIS 
arms) to match the lower resolution of the models in 
this regime. Finally, the observed spectrum is corrected 
for foreground Galactic redd ening (Ay = 0.0689) a dopt- 



ing the extinction maps of Schlegel et al. (|1998| and 



put in the rest-frame. The model spectra are convolved 
with a range of velocity dispersions that accounts for the 
measured velocity dispersion and associated uncertain- 
ties (particularly relevant in the optical range and for 
absorption index measures). 

In order to put constraints on the above mentioned pa- 
rameters, we combine information from individual stellar 
absorption features (see Figure 3]) with UV and optical 
colors. The indices we analyze here are chiefly sensitive 
to age and metallicity and minimally sensitive to dust 
and a/Fe abundance ratio (this is import ant because our 
models have solar-sc aled abundances, e.g. Worthc y|1994 



Thomas et al.|[2003 ) while the colors give constraints on 



dust attenuation. Specifically, we interpret the stellar ab 
sorption indices L>4000„, HS A + Hj A and [MgFe]', the 
FUV (1300 - 1780A) to NUV (1770 - 2730A) flux ratio 
and the observer-frame J-H color. The parameter esti- 
mates and associated uncertainties (median and 16 th and 
84 th percentiles of the PDF, respectively) are summa- 
rized in Table [6] We derive a luminosity- weighted (mass- 
weighted) age of 0.8 Gyr (1.2 Gyr) with an accuracy of 
<30%. The rather weak and low-S/N Mg and Fe features 
do not allow us to put strong constraints on stellar metal- 
licity which is estimated to be \og(Z*/Z & ) — 0.02+ ; 41 . 
The relatively large errors on stellar metallicity limit our 
ability to accurately estimate dust attenuation, however 
the rather strong 4000A-break and Balmer absorption 
lines in combination with the faint and flat UV require 
a fair amount of dust in this galaxy. The attenuation 
by dust is estimated to be Ay = 0.77i ; 32 in the optical 

and Apuv = 2.5lg'g in the UV. The attenuation in the 
optical reflects the optical depth in the ISM, HTy, experi- 
enced by stars older than 10 7 yr, while the attenuation in 



the UV traces well the total effective optical depth ry of 
young stars. Finally, the stellar mass, obtained by nor- 
malizing the models to the observed luminosity between 
the rest-frame 5000-5800A (accounting for dust attenu- 
ation), is estimated to be log(M*/M ) = 11.37±g;io- 

We tested how sensitive the derived parameters are to 
the observational constraints adopted and found good 
agreement among the various age estimates and on the 
need for a significant amount of dust. We note though 
that the inclusion of UV constraints is particularly im- 
portant for estimating dust attenuation. Our default fit 
combines the power of absorption features in alleviating 
parameter degeneracies and the sensitivity of UV and 
optical colors to dust attenuation. We take the standard 
deviation of the physical parameters obtained fitting dif- 
ferent observational quantities as an estimate of the pos- 
sible systematic uncertainties. The results are summa- 
rized in Table ?? while the comparison between different 
fitting methods is discussed in Appendix [A] 

The estimated mass- weighted age of 1.2 Gyr implies 
that the galaxy must have formed at z > 3 [zt = 3.2+Qg 
), while the luminosity- weighted age implies that star 
formation continued at least for other 0.4 Gyr (i.e., until 
z <~ 2.7). We look for the presence of weak nebular emis- 
sion lines after subtracting t he bestfit stellar cont inuum 



using the PLATEFIT code ( Tremonti et al. 
analyzing the residuals we find a 3cr upper 
Ha flux of 3.5 x 10 _17 erg/s/cm 2 



2004Q By 



cutt 



lmit to the 
Adopting the Kenni- 



( pL998 ) calibration between SFR and Ha luminosity 
(SFTlp^/yr] = 5.2xlO~ 42 L(Ha)[erg/s] for the Chabrier 
IMF), this translates into an upper limit to the current 
SFR of 5.5M0/yr. If we correct the Ha luminosity for 



5 This code finds the bestfit non-negative linear combination 
of SSPs and the bestfit dust attenuation that describe the stellar 
continuum in regions free of emission lines. The bestfit stellar 
continuum (and any smoothed residuals) is then subtracted from 
the original spectrum in order to obtain the 'pure' emission line 
spectrum. 
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TABLE 3 

Physical parameter estimates and associated 
uncertainties obtained by fitting stellar 
absorption indices, the rest-frame uv and 

optical color. the parameter estimate is the 
median of the associated pdf, while the 

statistical uncertainty is measured from the 

16 th AND 84 th PERCENTILES OF THE PDF. THE LAST 
COLUMN GIVES THE SYSTEMATIC UNCERTAINTY 
ESTIMATED BY FITTING DIFFERENT OBSERVATIONAL 
QUANTITIES (SEE APPENDIX |A| FOR. DETAILS). 



Parameter 



Estimate Systematic uncertainty 



log{t r /yr) 


8.90 +0 ;S° 9 


0.05 dex 


log(t™/V) 


9.081°,;™ 


0.03 dex 


log(Z*/Z Q ) 


02+ 20 


0.06 dex 


log(M„/M Q ) 


j, o 7 +0.13 


0.12 dex 


T V 


i 7s +0.88 


0.33 




o 64+ - 32 

u - o4 -0.28 


0.21 


Afuv 


9 co+0.89 
-0.81 


0.53 mag 


A V 


n 77+O.36 
u -' '-0.32 


0.27 mag 



the inferred dust attenuation in the optical, the SFR is 
< lO.3M /yr and the specific SFR is < 4.4 x 10 -11 j/r -1 . 

Finally, we note that the derived stellar metallicity is 
consistent within the uncertainties with the metallicity 
of eq ually massive z = 0.1 galaxies (see |Gallazzi et al.| 
2005| . The luminosity- weighted age is yoimger than the 
local population, as expected, but, if we assume that 
the galaxy would evolve passively since the redshift of 
observation until today, its present-day stellar age would 
be fully consistent with the local age-mass relation. 



5. SURFACE BRIGHTNESS PROFILE FITS 

We fit the 2D surface brightness profile of the galaxy 
in the three NIR UDS broad bands J, H and K , using 



galfit ( |Peng et aL]|2002| ). We assumed a |Sersic| J 1968 1 
model and allowed the Sersic parameter n, the effective 



radius r e , the axis ratio ab and the position angle PA 
to vary freely. The point spread functions (PSFs) were 
estimated from high S/N stacks of stars extracted from 
the UDS survey (as described in |Williams et aL 2010 ). 
In Figure [l] we show the K-band image of the galaxy 
and illustrate the fitting process. There is a foreground 
(zphot = 0.58) galaxy - 3" NW of UDS19627, which can 
potentially influence the fit. We perform three different 
fits to estimate the uncertainty it introduces. 1: Model 
the two galaxies simultaneously, 2: mask out the fore- 
ground galaxy and fit the main galaxy, 3: model and 
subtract the secondary galaxy (masking out the main 
galaxy), and fit the main galaxy in cleaned image. In 
Table [4] we list the best fitting parameters of the main 
galaxy for the three fits. The parameters are slightly sen- 
sitive to how the secondary galaxy is taken into account 
in the fit, but agree within the formal 3cr uncertainties. 

We adopt the best fitting parameters of Fit 1 in the 
remainder of the paper, and include the standard de- 
viation of the parameters in the three fits in the error 
estimates, in addition to the formal fitting errors from 
galfit. These parameters are printed in bold face in 
Table [4] In the same way we also derive the best fitting 
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Fig. 3. — X-Shooter spectrum (binned to 3A per bin for illus- 
tration purposes) zoomed over the region of the absorption in- 
dices used in the fit (the two bandpasses defining the 4000A-break 
and the central bandpass of the other indices are indicated by the 
hatched regions). The vertical dotted lines indicate the main ab- 
sorption features detected. The red spectrum is the bestfit model 
to the absorption indices, the UV flux ratio and the observer-frame 
J-H color. 

parameters and errors in the H and J bands. The Sersic 
parameters agree in the three wavebands. The circular- 
ized sizes derived in the H and J band are slightly larger, 
but more uncertain and agree within 2a. 

5.1. MIR emission 

The UDS was observed as part of the Spitzer UKIDSS 
Ultra Deep Survey (PI Dunlop). At 24/zm, the posi- 
tion UDS19627 is blended with a bright source coincid- 
ing with the nearby z = 0.58 galaxy (red cross in Figure 
[I]). We model and subtract this source with galfit using 
a bright nearby point source as PSF. The galaxy is not 
detected in the residual image. In an 7.5" aperture we 
derive a flux of 0.1 ± 20^Jy, where we have included the 
typical error in this size aperture quoted in the Spuds 
survey documentation^ The non detection rules out 
significant amounts of obscured star forma tion or AGN 
activit y. Using the method of described in Franx et al. 
( |2008| we derive a 2a upper limit on the star formation 
rate SFR < AOA^yr" 1 . 

6. STELLAR VELOCITY DISPERSION 

The presence of several strong stellar absorption fea- 
tures in our spectrum allows for high quality measure- 
ment of the stellar velocity dispersion. This determina- 
tion is based on the assumption that the observed line 
broadening is dominated by disordered motions of the 



6 see http://irsa.ipac.caltech.edu/data/SPITZER/SpUDS/documentation/RE 
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TABLE 4 

Parameters of the 2D surface brightness distribution derived with galfit. a c is the 

EFFECTIVE SEMI-MAJOR AXIS. N IS THE SERSIC PARAMETER, b/A IS THE AXIS RATIO, 
r e>c = a 0ima j or * \fbja IS THE CIRCULARIZED EFFECTIVE RADIUS. FlT 1-3 IN THE THREE FIRST 

rows are done in the k-band, dealing with the foreground neighbor galaxy in 
different ways to estimate its effect on the derived parameters of the main galaxy. 
In Fit 1 the two galaxies are modelled simultaneously. In Fit 2 the main galaxy is 

FIT WITH THE NEIGHBOR GALAXY MASKED OUT. In FlT 3 THE MAIN GALAXY IS MODELLED IN AN 

IMAGE WHERE THE NEIGHBORING GALAXY HAS FIRST BEEN MODELLED AND SUBTRACTED. THE 
ERRORS QUOTED FOR THE THREE FITS ARE THOSE RETURNED BY galfit. THE NEXT THREE ROWS, 
LABELED "BEST FITS" , ARE THE PARAMETERS OF THE MAIN GALAXY FROM SIMULTANEOUS FITS IN 
THE K, H AND J BANDS. THE ERRORS QUOTED HERE INCLUDE THE STANDARD DEVIATIONS FROM 
OF THE BEST FITTING PARAMETERS FROM FlT 1-3. 



Fit Band 


Mag 


a e ["] 


n 


b/a 


x 2 


r e , c [kpc] 


1 K 

2 K 

3 K 


20.24 ±0.01 
20.18 ±0.01 
20.26 ±0.01 


0.44 ±0.01 
0.47 ±0.01 
0.44 ±0.01 


1.47 ±0.09 
1.67 ±0.07 
1.42 ±0.07 


0.74 ±0.02 
0.77 ±0.02 
0.73 ±0.01 


0.948 
0.963 
1.073 


2.77 ±0.06 
3.04 ±0.07 
2.70 ±0.05 


Bestfit with systematic errors 


1 K 

1 H 
1 J 


20.24 ± 0.04 

20.54 ±0.11 
20.99 ±0.07 


0.44±0.02 

0.48±0.05 
0.50±0.03 


1.47±0.15 

1.43±0.35 
1.46±0.23 


0.73 ±0.03 

0.81 ±0.03 
0.80 ±0.05 


0.948 

0.938 
0.937 


2.77 ±0.11 

3.24 ±0.37 
3.20 ±0.20 



stars and that the contribution due to bulk flows or ro- 
tation is minimal. We have therefore fit the spectrum 
using a set of resolution-matched stellar template spec- 
tra to model the effect of the velocity broadening. 

6.1. Penalized Pixel Fitting 

The p ublic penalized pixel fitting co de, PPXF, devel- 
oped by Cappellari & Emsellem ( 2004 ) , has been used in 
the literature to fit velocity dispersi ons of high redshift 
galaxies (e.g.,|C appclla rTet al.|2009||Onodera et al.|2010 
van de Sande et al.||201ll ). We have run pPXF with the 
MILES empirical stellar library ( Sanchez-Blazquez et al. 



2006[|Falc6n-Barroso et al.|2011[ ). As input, pPXF takes 
the raw X-Shooter spectrum (excluding the observed K- 
band) along with all 985 stellar spectra from MILES. 
We have shifted the observed spectrum to the rest-frame, 
smoothed it with a Gaussian to match the instrumental 
resolution of the templates and used logarithmic binning 
to put all spectra on the same wavelength grid. The 
code then constructs an optimized linear combination of 
a small subset of the 985 spectra, in this case, nine tem- 
plates are chosen and combined. Using only this sub- 
set we have run 1000 fits where the target spectrum is 
randomized within the errors and the weights for the 
nine input templates are allowed to vary (including zero 
weight). The resulting distribution of stellar velocity dis- 
persions are shown in Figure [4j The best fitting velocity 
dispersion from this analysis is er* = 318 ± 53 km s . 

The inferred velocity dispersion changes in relation to 
the radii at which it is measured. Using Equation 1 from 
Cappellari et al. (2006) and assuming that the scaling 
from the spectrum to the broad-band fluxes is only due 
to missing light from the profile, we infer the correction 
to er at r e to be only a few percent. 

We performed an independent test of the derived ve- 
locity dispersion through a chi-squared analysis in which 
we fitted gaussian convolved, stellar template spectra of 
five different types (AV, F5V, G8IV, K0III, K0V) ob- 
served with X-shooters UVB arm, which corresponds ap- 
proximately to the same restframe wavelength range as 
sampled by the NIR arm for the z=2 galaxy. The AV 




<a t > =318 ±53 km s~ 



MlJ 



300 350 

a* [km s _1 ] 

Fig. 4. — Histogram of best-fit velocity dispersions for the 1000 
randomized realizations of the pPXF fitting. Each fit is produced 
by a new linear combination of the initially chosen nine stellar 
templates. 



and F5V stars and linear combinations hereof provide 
good fits with velocity dispersions in the range 260 — 300 
km s _1 with typical uncertainty of 50km s _1 , in good 
agreement with the results from pPXF. The analysis is 
described in detail in Appendix [B| 

6.2. Dynamical mass 

We estimate the dynamical mass of the galaxy from 
the derived velocity dispersion and effective radius (in 
the if-band): 

M dyn = f3r e a 2 JG, (1) 

where /3 is a constant that depends mainly on the struc- 
ture of the galaxy. A value of = 5 is co mmonly used in 
the literature (e.g. J0rgensen et al.|1996 ). From detailed 
analysis and modeling ot spatially resolved kinematic ob- 
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servations of a sample of 25 local E and SO galaxie s and 
comparison to det ailed modeling of stellar orbits, |Cap- 
pellari et al. (2006) found that (3 — 5.0±0.1 accurately re- 
produces galaxy dynamical masses. Adopting this value, 
we find M dyn = 3.26 ±0.99 x lO n M . From simulations 
it is predicted that /3 will depend on the struct ure of 
the galaxy, increasi ng with smaller Sersic n values ( Cap- 
pellari et al. 2006). For the best fitting n = 1.47, a 
value [3 ~ 7.5 is predicted. If we adopt this value we 
find M dyn = 4.88 ± 1.49 x lO n M , but as this depen- 
dency on structure has not been verified observationally 
( Cappellari et al.]|2006 ), we will adopt (3 = 5 in this pa- 
per, which also makes compariso n to the literature more 



straight forward (but see Section 7. 



6.3. Spectral Classification of UDS19627 

The spectrum of UDS19627 has the characteristic 
"shark-tooth" shape of a post starburst galaxy, consis- 
tent with the low derived SFR, and constraints on the 
stellar age and the prominent Balmer series absorption 
lines. The derived magnitudes (see Table El) corresponds 
to an Lk ~ ga laxy (assuming 7l* =2 = 21.3, 

Marchesini et al. 20071. The Sersic fit to the surface 
brightness distribution yields a low n = 1.47, also more 
similar to what is found in late type (and post starburst) 
galaxies, than what is found in local early types. 



7. 



SAMPLE OF Z ~ 2 GALAXIES WITH MEASURED 
VELOCITY DISPERSIONS 



Prior to the results presented here, velocity dispersions 
of three quiescent, massive z ~ 2 galaxies have been pub- 
lished. We include the properties of these galaxies in 
our analysis. Properties of their stellar populations have 
all been derived assuming a Chabricr IMF. To minimize 
systematics errors, we recalculate the dynamical masses 
using effective radii derived from restframe optical, and 
assuming the same value for /3 = 5. 

7.1. COSMOS-254025 

Onodera et al. ( 2010[) reports cr* < 326 km s _1 derived 
from a Subaru/ MU1RCS spectrum of a z — 1.82 galaxy 
selected from the COSMOS survey. This is an upper 
limit due to the limited resolution of the spect rograph. 
Using the re sults of surface brightness fits by IMancini] 



et al. 



(20101 (r e = 5.79 ± 0.61 kpc and n = 4.14) from" 
J5F 



HST/ACS F814W band observations, they derive a dy- 
namical mass of Al < 7 x 1O 11 M . For this analysis we 
measure the effective radius in NIR images from the Ul- 
travista survey (McCracken et al, in prep.), sampling the 
restframe optical, which is where the velocity dispersion 
is measured, and is a better tracer of the stellar mass 
than the ACS I band (restframe UV), as z ~ 2 galaxies 
can have strong morphological k-corrections (e.g. |Toft et| 
al. 2005 Cameron et al. 2011). Indeed, comparison ot 
the ACS and Ultravista images shows that the galaxy is 
extended and clumpy in the restframe UV, but smoother 
and more centrally concentrated in the restframe optical. 

Using an identical approach as described in Section [5] 
we run galfit on the deep J, H, and K band UltraVISTA 
survey observations (McCracken et al, in prep.), mask- 
ing out neighboring galaxies, and using a bright nearby 
star as PSF model, we find a mean effective radius of 
r e = 2.4 ± 0.4 kpc, and n — 2.5 ± 0.4. The derived 



et al. 



60% smaller than the value derived by Mancini 
(|2010|) in the F8UW band. While the galaxy is 



only margmaly resolved in the UltraVISTA data, sizes 
measured for galaxies using the procedure followed have 
been demonstrated to be reliable for galaxies with intrin- 
sic sizes several times smaller than the PSF (|Trujillo"et] 



al.| |2006| Toft et al.p009| ) when the galaxy is detected 
at high S/ IN and a good PSF model is available, as is 
the case here. Note that the extend of the galaxy in 
the ACS image is twice the size of the Ultravista PSF 
(FWHM~ 0.7"), so if the restframe optical morphology 
had been similar it would have been well resolved in the 
Ultravista images. 

With our measured r e the upper limit on the dynami- 
cal mass becomes Md yn < 2.97 x 10 11 M Q . From the sur- 
face brightness fits to the UltraVISTA data we derive an 
H-band magnitude of 18.42 ±0.01 and K = 19.81±0.01, 



cor respon ding to an Lk = 22L* K galaxy (Marchesini et 
aL]|2007|). From a com bined fit of the spect rum and 



Broad band photometry Onodera et al. (2010) quote a 



stellar mass of 3-4xlO n M0, and a luminosity weighted 
age of 1-2 Gyr, corresponding to a formation redshift of 
Zf = 2.5 — 4. The shape of the spectrum is typical of a 
post starburst galaxy, with prominent Balmer absorption 
features. 

7.2. NMBS-C7447 
From an X-shooter spectr um of a redshift z = 1.800 



galaxy in the COSMOS field, van de Sande et al. (20111 



294 ± 51km s" 



(1.7 ± 0.5) x lO n M , and 



derive a velocity dispersion of a 
a dynamical mass M ( i yTl 

from HST/WFC3 observations in the F160W band, an 
effective radius of r e = 1.64 ± 0.15 kpc and Sersic pa- 
rameter n = 5.3 ± 0.4. The X-shooter spectrum has 
a typical post starburst shape, with prominent Balmer 
lines, and the stellar population synthesis fits to the 
spectrum is consistent with little ongoing star formation 
(O.OO2Af yr" 1 ), a relatively young age 0.4 Gyr, corre- 
sponding to a formation redshift of Zf ~ 2 and a stellar 
mass 1.5 x 10 n M Q . In the MUSYC catalog ( |Gawiser et] 
|al.|2006[ ) it has H = 19.85 and K = 19.75, corresp onding 
i Lk ~ 23L^ galaxy ( |Marchesini et al.||2007| ). 



to an Lk 



7.3. MUSYC 1252-0 



van Dokkum et al. 



(2009) report a a = 510Z g5 kms 



derived from a deep Gemini/ GNIRS spectrum of a z = 
2.186 galax y selected from the MUSYC survey ( |Gawiser 
et al.||2006|). It has H = 21.31, K = 21.03, corr espond- 



mg to a L=7L K galaxy (Marchesini et al. 2007). Based 
on surface brightness fits to HST/1N1C2 1160VV observa- 
tions, they find r e = 0.78 ± 0.17kpc and n = 3.2 ± 0.9 
van Dokkum et al. 2008). From these numbers, and 



assuming /3=5, we calculate M, 



dyn 



= 2.4 ±^_x 10 n tf B 



whi ch is i dentical to the value quoted in |van Dokkum et| 
al. (2009). The shape of the GNIRS spectrum is consis- 
tent with old, evolved, quiescent stellar population, with 
little ongoing star formation (SFR^ 1 — SMQ-yr^ 1 ) and 
age 1.3-2.0 Gyr, corresponding to a formation redshift 



of Zf ~ 4- 7 ( Kriek et al.|2009 |, quite different from the 
younger, post starburst spectra of the other three z ~ 2 
galaxies considered. This galaxy also has a significantly 
larger velocity dispersion and smaller size than the other 
three galaxies (see Figure [5|. 
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TABLE 5 

Properties of the sample of z ~ 2 quiescent galaxies with measured velocity dispersions. 



ID 


z 


re 
[kpc] 


n 


log(M») 
[M Q ] 


log(Agc) 
[Gyr] 


<*obs 
[km s- 1 ] 


log(M dyn ) 
[Mq] 


[km s- 1 ] 


UDS 19627 
NMBS-C7447 


2.0389 
1.800 


9 77 +0.11 
■ 2 -' '-0.11 
1 64+ - 15 


1 47+0.15 
LAI -0.15 

r on+0.40 

o.ou_ 40 


11 oy+0.13 
- L1 - o '-0.10 
ii ,0+0.20 
■ U - L °-0.20 


a nn+o-io 
8-90_ 09 

8.60 


318l|| 
294±|J 


11 51 + 011 
iJ -- oi -0.11 

11 23 + - 11 


339 
281 


COSMOS-254025 


1.820 


2.40t°;f 




11 54+006 


Q 1B+030 


<326 


<11.47 


354 


MUSYC 1252-0 


2.186 


0.78lg;i? 


o 9n +0.90 
°' zu -0.90 


11 36+ 04 


9 24+ 10 
3 - z -0.13 


510t 9 f 


1 1 07+0.25 
1J -' 3 '-0.13 


469 



7.4. Comparison to Low Redshift Galaxies 
In Figure [5] we compare the properties of the four z ~ 2 



galaxies with measured stellar velocity dispersions (col- 
ored points) to local z < 0.2 early type (black contours), 
and post starburst (cyan contours) galaxies from the 
SDSS New Yo rk University Value-A dded Galaxy Catalog 
(NYU-VAGC |Blanton et al.||2005[ ). The former galaxies 
are selected using the criteria: SFR< 1 M^yr -1 , and 



Sersic n > 3, the latter using the criteria of Goto (2005 1 . 
Effective radii are measured in the z-band. 

Three of the z ~ 2 galaxies (UDS19627, NMBS-C7447 
and COSMOS-254025, represented by stars) follow a 
mass-size relation with similar slope as for local early 
type galaxies, but offset to smaller sizes at a given mass 
(in the following we will refer to these three galaxies as 
the post starburst SEEDs). One galaxy, MUSYC 1252-0 
(square), is an outlier with respect to these (we refer to 
this galaxy as the evolved SEED). 

In Figure [5] (a) we show the relation between r e and 
stellar mass. The z ~ 2 post starburst galaxies follow 
a relation similar to the local early type relation (solid 
line), but offset by a factor of 3.3 ± 0.2 to smaller r e for 
a given stellar mass (dashed line). This offset is similar 
to what has previously been found for larger samples of 
stellar mass selected quiescent galaxie s at z ~ 2 with 
photometric masses and redshifts, (e.g. Toft et al.|2009 1, 
but smaller than the factor of five offs et found for color 
selected z ~ 2 quiescent galaxies (e.g. Toft et aL"1|2007[ 
Zirm et al][2007l [Cimatti et al.|[2008l Ivan Dokkum et| 
^2008), exemplified in these plots by the color selected 
M U S Y 1252-0 galaxy. The stellar mass-size relation of 
low-z post starburst galaxies is on average shifted along 
the early-type relation to larger masses and sizes. In 
Figure (b) we show the same relation for dynamical 
mass versus effective radius. The z = 2 post starburst 
galaxies follow a relation similar to the stellar mass-size 
relation, but with a smaller offset and scatter. r e for the 
z ~ 2 post starburst galaxies are on average 2.5 ± 0.2 
times smaller than local early type galaxies of similar 
dynamical mass. As illustrated in Fig[5](c) the difference 
in offset in (a) and (b) is due to a difference in the derived 
stellar to total mass ratio. 

For local early type galaxies there is a tight correlation 
between Md yn and M„, with the ratio of dynamical to 
stellar mass (within r e ) increasing fro m the least massive 
to the m ost massive gala xies ( see also [Pa dmanabhan et 
aT| |2004| |Gallazzi et aL"1|2006| |Taylor et~al.||2010p . UrT 



der the assumption of structural homology, this can be 
interpreted as an increasing dark matter fraction with 
mass. The local galaxies with \og(Mdy n /M Q ) ~ 10.5 
are thus baryon dominated, while the most massive 



log ( Mdyn /Mq) > 11.5 have much larger dark matter 
fractions. A similar trend has been shown to be valid 
in early type galaxies out to z 



1 (e.g. 



2005) |Rettura et al |[2006| ) . While the ,: 
consistent with the local re 



Ferreras et al. 
2 galaxies are 
relation within the error bars, it 
is striking that the best fitting values all fall close to the 
dashed line (M* = Mdyn) indicating a lower dark matter 
fraction than in similar mass early type galaxies at lower 
redshift. We explore this further in Figur^5](d) where we 
plot 1 — M* /Mdyn'- a proxy of dark matter fraction as a 
function of dynamical mass (assuming homology). The 
central dark matter fraction of the z ~ 2 galaxies scatter 
around a mean value 1 — M*/Mdyn = 0.18 ± 0.20, while 
the low redshift galaxies in the same dynamical mass 
range have a mean of 0.46 ± 0.23. Note from (b) that the 
dynamical mass-size relation of the local post starburst 
galaxies is shifted from the early type relation to larger 
sizes at a given mass. From plot (c) and (d) it can be 
seen that the low dark-matter fraction found in the z ~ 2 
galaxies is consistent with what is found in local post 
starburst galaxies, which show a large range in M*/Mdym 
but on average have Md yn ~ M* . The observed low cen- 
tral dark matter fraction of the z ~ 2 galaxies compared 
to that in local early type galaxies could have implica- 
tions for the evolutionary path to lower redshift. We 
explore this further in the discussion (Section [8|. 

In Figure [5] (e) we plot the effective radius versus the 
velocity dispersion. The average velocity dispersion of 
the z ~ 2 post starburst galaxies is 1.8 ±0.5 times larger 
than in low redshift galaxies of similar size, while the 
velocity dispersion of MUSYC 1252-0 is 4.1 ± 2.1 times 
larger than local galaxies of similar size. The local post 
starburst galaxies have a similar distribution as local 
early types but shifted to larger sizes and velocity dis- 
persions. 

In Figure [5] (f) we plot the inferred (dynamical) mass 
density within the effective radius S = 0.5Md y n/(^rl), 
versus the dynamical mass Md yn . The z ~ 2 post 
starburst galaxies have average surface mass densities 
5.8 ± 1.2 times higher than the mean in local early types 
of similar mass. The relation for the local post starburst 
galaxies are shifted to smaller surface mass densities. As 
in the dynamical mass-size plot, the offset of the z=2 
post starburst galaxies from the local post starburst re- 
lation is similar to the offset of the evolved z=2 galaxy 
from the local early type relation. 

7.5. Evolution Through Merging? 

Dry merging has been suggested as a mechanism to 
"puff up" the SEEDs and evolve them into agreement 
with the mass-size relation observed at low redshift 
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Fig. 5. — Properties of the four z ~ 2 massive quiescent galaxies with measured velocity dispersions, compared to early type and post 
starburst galaxies in the local universe | |Blanton et al.|2005| drawn from the SDSS NYU-VAGC). The blue star is the galaxy considered in 
this paper. The two red stars and the red square are tne post starburst and evolved SEEDs from the literature. The black/cyan contours 
are low redshift early type/post starburst SDSS galaxies. The full black lines are linear fits to the local galaxies, the dashed lines are linear 
fits to the high redshift points, with the slop e fixed to th e local values. The orange arrows represent the evolution due to merging predicted 
derived from the cosmological s imulations of |Oser et al.||2012fl . The colored lines show the effect on the results of usin g the expression for 
P{n) of Cappcllari et al. (2006) when calculating the dynamical masses rather than assuming /3 = 5 (see Section |7.8| l. Top left: Stellar 
mass vs effective radius lop right : Dynamical) mass vs effective radius. Middle left: Dynamical mass vs stellar mass. Middle right 
Central dark matter fraction (1 — M* /M^yn) vs dynamical mass. Bottom left: velocity dispersion vs effective radius, Bottom right: 
Surface (dynamical) mass density vs dynamical mass. 
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2009 Newman 


et al.||^01l|. 


Previous studies oi the stellar mass-size re- 



lation have not been able to unambiguously distinguish 
between different merging scenarios, but minor merging 
is the most promising process. With the added dynami- 
cal information we can now start to test different merger 
scenarios in greater detail than what is possible with just 
stellar masses and sizes. The low inferred dark matter 
fraction of the z ~ 2 galaxies, compared to early types 
at low redshift, supports the merging scenario, as merg- 
ing can redistribute the dark matter within the effectiv e 



radius (Boylan-Kolchin et al. 2005 Oser et al. 2012) 



extended than the stars, increasing r e by adding stars 
in the outskirts can lead to an increased measured dark 
matter fraction by up to a factor of 2-3 because areas 
with higher dark matter content is included within the 
effective radius (Hiltz et al, in prep., T. Naab, private 
communication) . 
Base d o n simple virial a rguments [Bezanson et aL| 



Furthermore, if the dark matter profile is more radially 



(120091) and INaab et al.|(|2009|) show thatlny merging can 
lead to an increase in size of the remnant, and that mi- 
nor merging is likely to be the dominant process, as this 
is the most effecient process for size growth, in terms of 
added mass needed for the observed size evolution (about 
a factor of 2). Cosmological simulations studying the ef- 
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feet of merging on the size evolution of massive galaxies, 
show that dominant mode of mass growth is minor merg- 



ing, w ith a typical mass weighted ratio of 1:5 (Oser et al. 
2012|. The average integrated size growth between z=2 
and for M ~ 10 11 Mq galaxies in these simulations is 
a factor of 5-6, and an average mass growth of a factor 
of 2.1. The velocity dispersion in the remnants are ap- 
proximately 30% lower than in their progenitors at z=2. 
This evolution is almost identical to what is predicted 
from the simple virial approximations for a similar mass 
growth. Interestingly, in these simulations galaxies at 
z=2 typically have low dark matter fractions at z = 2 
(0.1-0.3, similar to what is derived here for the observed 
z ~ 2 galaxies), and these increase with time as a conse- 
quence of the merging by approximately a factor of 2. In 
Figure [5] we plot with a rrows the evolutio n from the cos- 
mological simulations of |Oser et al.| ( |2012[ ). For the z ~ 2 
post starburst galaxies, these predictions provide a good 
match to the observed evolution of the stellar mass and 
velocity-size relations. However, they fail to reproduce 
well the evolution in the dynamical mass-size and dy- 
namical mass density - mass relations where the galaxies 
end up above and below the local relations, respectively, 
as a consequence of the evolution of dark matter frac- 
tion in the simulations not being sufficient to account for 
the observed evolution. For the evolved z ~ 2 galaxy 
additional structural evolution is needed. Some of the 
additional evolution needed in the dark matter fraction 
may be attributed to merging with galaxies with higher 
dark matter fractions than those included in the simula- 
tions. In the local universe dwarf galaxies can be heavily 
dark matter dominated, with stellar masses ~ 10 7 M Q , 
and dynamical masses ~ 10 8 — 1O 9 M , corr esponding to 
dark matter fractions as high as 95% (e.g. Carignan fc l 
Beaulieu|[l989"l [Persic et aL||1996| |Strigari et"lil.j|2010[) . 
If the z ~ 2 galaxies grow via minor merging with pri- 
marily dark matter dominated dwarfs, it could be an ef- 
ficient way of increasing the dark matter fraction in the 
remnant. 

7.6. The Fundamental plane at z ~ 2 

From the measured velocity dispersions, effective radii, 
and the mean restframe surface brightness within the 
effective radius 

(AO = K-51og(r e ) + 2.51og(2^)-2.51og((l + z) 4 )-A v , 

(2) 

we can for the first time probe the fundamental plane 
of z ~ 2 quiescent galaxies: 



log(r e ) = alog(a) - /31og(I e 



(3) 



where V z is the ( AB) magnitude in the redshift ed V band, 
which following |van Dokkum fc Franx ( 1996J ) can be es- 
timated simply as V z = H + 2.51og(l + z) (since the ob- 
served ii-band resembles closely the 1^-band redshifted 
to z ~ 2), r e is the (TJ-band) effective radius, a is the 
velocity dispersion in km s _1 and I e — lO^ 4 *^). 
Here we adopt a = 1.2, and ft = 0.83, as derived by 



J0rgensen et al. (1996) for a sample of 225 galaxies in 
10 nearby clusters. Note that this quite strong assump- 
tion for the 3D orientation of the fundamental essentially 
builds in homologous evolution in the following analysis, 
but its the best we can do given the paucity of z ~ 2 
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Fig. 6. — Comparison of the restframe V-band fundamental plane 
of the 2 ~ 2 galaxies compa red to in the z=0.02 Coma cluster (black 
squares. iLuce y et al. | 1991b. and in the z=0.33 cluster CL1358+62 
(black circles, Kelson ct al. 2000). The stars are the post starburst 
SEEDs, the blue is tne one considered in this paper, the square is 
the evolved SEED. The z = 2 galaxies are offset significantly from 
the local relation. The lines indicate the effect of luminosity and 
structural evolution considered between z=2 and 0. The dotted and 
dashed lines are the evolutio n pr edicted by the minor and major 
merging scenarios in Section |7.4| and the full line is the effect of 
passive evolution of the M/L predicted by the best fitting stellar 
population synthesis models between their observed redshift, and 
today, assuming the best fitting formation redshift Zt, a Chabricr 
fMF and solar metallicity. These combined evolutionary effect are 
able to bring the evolved SEED close the local fundamental plane, 
but cannot fully account of the offset of the post starburst SEEDs 



data points. 

In Figure [6] we compare the rest-frame V-band fun- 
damental plane of the z ~ 2 galaxies to the funda- 
mental plane observed loc ally in the z ~ 0.02 Coma 



cluster (Lucey e t al 



ter CL1358+62 (Kelson et al. 2000) 



1991[) and in the z ~ 0.33 clus- 
The galaxies have 

been corrected for cosmological (1 + z) 4 surface bright- 
ness dimming. The slope of the z — 0.33 fundamental 
plane is similar to in the Coma cluster, but slightly offset. 

While it is not possible to derive robust estimates of 
slope or scatter of the z ~ 2 galaxies, we note that they 
show a significantly larger scatter than the local relation, 
and a large offset from it. As in Figure [5] the evolved and 
the post starburst z ~ 2 galaxies fall in different parts of 
the plot. Assuming the local slope, the post star burst 
galaxies are on average offset by 1.6 ± 0.2 dex from the 
local relation. The evolved SEED is closer to the local 
relation. 

As the velocity dispersion traces the total mass, and 
the surface brightness traces the light, the offset between 
the fundamental planes at different redshifts is often in- 
terpreted as an offset in mass-to-light ratio caused by 
passive aging of the stellar populations, which shifts the 
plane up, as the surface brightness fades, while the ve 
locity dispersion stays the same (e.g. van Dokkum & van 
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der Marel|[2007| . 

'1'he offset observed for the z = 2 post starburst galax- 
ies is too large to be explained by passive evolution. The 
solid lines represents the maximum evolution of the sur- 
face brightness that can be attributed to passive aging of 
the stellar population between z=2 and z=0, assuming 
solar metallicity, a Chabrier IMF, and that the stars were 
formed in a single stellar burst at their derived Zf. We 
note that this relies on the ages being correct. If these 
are systematically overestimated there could be room for 
more passive evolution. It is clear that in order to end up 
on the local relation the galaxies also have to go through 
structural evolution. The minor merging scenario de- 
scribed above would also have an effect on the funda- 
mental plane as it changes the velocity dispersion and 
size, and thus also the surface brightness of the galaxies. 
The expected shift of the galaxies due to this process is 
indicated with dotted lines. 

These evolutionary effects in concert brings the z=2 
post starburst galaxies closer to the local relation but 
still systematically under. The offsets range from 0.2-1 
Dex (with a mean of 0.5±0.4). Interestingly, the evolved 
SEED ends up on the local relation through these pro- 
cesses, however from Figure [5] we know that it needs to 
go through additional structural evolution to end up on 
the local mass-size relation. 

A possible origin of the additional offset from the local 
fundamental plane that can not be accounted for by pas- 
sive evolution and minor merging is that the evolution 
between z=2 and may be non homologous. This is not 
unlikely as we know that at least some of the z = 2 galax- 
ies have to change their surface brightness profiles from 
exponential disk to deVaucouleurs profiles. Local post 
starburst galaxies have in some cases been shown to be 
offset by a similar amoun t from the funda mental plane 
of local elliptical galaxies ( Yang et al.|2004 ), which is at- 
tributed to the very young stellar ages of those particular 
galaxies. At the older observed ages of the z ~ 2 galaxies 
considered here, this is however not a very likely expla- 
nation. Metallicity and/or age gradients in the z ~ 2 
galaxies could be a potential source of uncertainty, how- 
ever higher resolution multicolor imaging is needed to 
adress this. 

7.7. Clues to the formation mechanism of SEEDs 

The observations presented here hold important clues 
to the formation scenario of SEEDs which is still not 
well understood. The uniformly old and extremely com- 
pact stellar populations suggest that the majority of the 
stars formed in a major nuclear starburst at high red- 
shift. Simulations indicate that highly dissipational pro- 
cesses on short timescales are plausible mec hanisms for 
creat i ng the compact stellar populations (e.g. Naab et al. 
2007 2009). A poss ible scenario is ma jor gas rich merg- 
ers at high redshift ( Wuyts et al.|2010" ), in which the gas 



the density of stars in SEEDs |Greve et aL 2005 Tac 



coni et al.||2006[ |2008| |MichaIowski et al.|[20TTj T The 



is driven to the center, igniting a massive nuclear star- 
burst (> lOOOA/Qyr- 1 ), followed by an AGN/QSO phase 
which quenches th e star formation, and leaves behind a 
compact remnant (Ho pkins et al.|2006 ). Sub-mm galax- 
ies (SMGs)/Ultra lummous infrared galaxies (ULIRGS) 
may be examples of this process in action. Observations 
show that most of them are major mergers, with large 
amounts of dust enshrouded star formation, and high 
central concentrations of molecular gas, comparable to 



compact structure, post starburst nature and (at least 
in the case of the galaxy studied here), relatively high 
extinction, make them very likely descendants of z > 3 
dusty SMGs. With the accurate ages obtainable from 
NIR spectroscopy of quiescent z ~ 2 galaxies, we can 
begin to constrain the number densities of their progeni- 
tors as a function of redshift, e.g. if one assumes a direct 
evolutionary link between quiescent z ~ 2 galaxies and 
z > 4 SMGs their relative number densities can be used 
to co nstrain the duty cy cle of starformation to be ~ 50 
Myr ( |Capak et al^|2008f 



7.8. Caveats 

As the z ~ 2 galaxies are spatially unresolved in the 
spectroscopic observations, it is not possible to determine 
if some of the broadening of the absorption lines, which 
we interpret as velocity dispersion, could be due to ro- 
tation. This may be a likely scenario given that some of 
their surface brightness profiles are best fitted by expo- 
nential disk like profiles, and there is evidence from recent 
high resolution observations that a significant fractio n of 
SEEDs may have fla ttened disk like morphologies 



van 



der Wei et al. 2011). This could lead to a systematic 
bias in their velocity dispersions, and dynamical masses. 
This, in turn, would mean that the derived offset from 
the local dynamical mass-size relation could be biased. 
Other sources of possible uncertainties in the dynami- 
cal masses are that the assumptions of homology and 
isotropy, needed to calculate the dynamical mass from 
the velocity dispersion may not be fully valid. Some stud- 
ies have suggested that the sizes of compact z ~ 2 galax- 
ies could be underestimated due to faint profile wings 
not detect ed due to cosmologic al surface brightness dim- 
ming (e.g. |Mancini et al.|2010 l, which would lead to the 
dynamical masses being under estimated. However, con- 
sidering the bright magnitudes of the galaxies considered 
in this paper, and that ultra deep studies tracing the sur- 
face brightness profile of compact z ~ 2 gala xies out to 
> 10r e have failed to detect such wings (e.g. |Szomoru et| 
aT||2010| 2011), this is not likely to be a very significant 
effect compared the other potential systematic errors. 

The derived low dark matter fraction is sensitive to the 
(3 = 5 assumption made when calculating the dynamical 
masses. From theory it is predicted that j3 should de- 
pend on the structure of the galaxies in a way that can 
be approximated through the following dependency on 
the sersic n parameter : f3 = 8.87 - 0.831n + 0.024n 2 



( Cappellari et al 
shows the effect o 



2006). The colored lines in Figure [5] 
assuming this P(n) correlation. The 
dynamical masses increase 3-50% for the z ~ 2 galaxies, 
resulting in a higher me an average dark matter f ractions 
(0.29 ± 0.18). However, [Cappellari et al.| d2006| ) showed 
that there is no observational support tor this correlation 
and real galaxies (with n=2-10) are best fitted by f3 ~ 5, 
we therefore adopt this value. 

Another potential uncertainty is the Initial Mass Func- 
tion (IMF) assumed in the stellar population synthesis 
modeling. Theoretical and observational studies have ar- 
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and z = galaxies may introduce a bias in the esti- 
mated stellar masses. Observations of z < 1 massive 
early type galaxies acting as gravitational lenses sug- 
gest that a Safpet er IMF provides a better fit than a 
Chabrier IMF (e .g. |Grillo et~aT1[2009l | Auger et al.plO 
2011). Assuming a Salpcter IMF rathci 



Sonnenfeld et ai.lMll). 
j-i — „„ / 'k„u ■„! i \ \ 6 



salpeter IMF rather 
than an Chabrier IMF for the z ~ 2 galaxies would lead 
to approximately 1. 75 times higher stellar masses (e.g. 
Gallazzi et al. 2008), resulting in stellar masses signifi- 
cantly higher than the dynamical masses, strongly dis- 
favoring a Salpeter IMF for the z ~ 2 galaxies. The 
derived stellar masses of the z ~ 2 galaxies has a rela- 
tively weak dependance on whether or not the IMF is 
bottom light, at the ob served stellar ages (~ 1 — 2 Gyr, 
Marchesini et al. 2009), but at the older ages of local 
galaxies (~ 10 Gyr) the derived stellar mass is more sen- 
sitive to the underlying IMF, leading to a potential bias 
between the z = 2 and z — mass size relations. As a 
consistency check, we calculate the so called "inferred" 
velocity dispersion, from the stellar mass and effective 
radius: a? n f — GM*/f3r e . The derived a in f of the z ~ 2 
agrees well (within ~ 10%) with the measured velocity 
dispersions (as expected since M dyn ~ M* , see Table [5f. 

As evident from Figure [T] UDS19627 is located close 
to a foreground galaxy (angular separation 9 = 2.83"). 
This galaxy is responsible for a gravitational lensing ef- 
fect that causes a bias in the derived brightness, derived 
stellar mass, size, and dynamical mass of the z ~ 2 back- 
ground galaxy. To estimate this effect we construct a 
simple lensing model of the system. We assume for the 
foreground galaxy a photometric redshift zi — 0.58 and 
stellar mass log(M* iZ /M Q ) = 10.32±3;o5, where both val- 
ues are derived from broadband B, r ,i ,z, J and K band 
photometry, using the same model library as for our main 
target. We then compare this stellar mass value to those 



of th e SLACS lens galaxies (Bolton et al. 2006 Grillo 



20101 to estimate the effective velocity dispersion a sis 
of a singular isothermal sphere (SIS) model that we adopt 
here to describe the total mass distribution of the lens. 
We find that a value of 145 km s" 1 would be typical for 
lens galaxies with stellar mass values similar to that of 
our foreground galaxy. For a SIS model, we recall that 
the Einstein radius 0Em is equal to 



Ei n = 47r(cr S7S /c) 2 D ls /D 



(4) 



where c is the speed of light and Di s and D os are the an- 
gular diameter distances between the lens and the source 
and the observer and the source, respectively. The mag- 
nification factor /i at an angular distance 9, that is larger 
than 9Ein, from the center of the lens is given by 



n{9) = l + 9 Em /{9-9 Em ). 



(5) 



Following the previous two equations, we calculate a 
magnification factor of 10-20% for the z ~ 2 galaxy. The 
brightness and stellar mass are thus over estimated by 
this amount and the size by the square root of this. By 
applying this correction, we would get a 0.07 dex lower 
stellar mass and a 0.03 dex lower dynamical mass. Given 
the relatively small weight of the lensing effect compared 
to the other uncertainties in our analysis and the approx- 
imations adopted in the lensing modeling, we decide not 
to correct our results for the lensing effect. We note that 
lensing is likely to be a general issue that needs to be 



taken into account when analyzing samples of the very 
brightest z ~ 2 quiescent field galaxies. 

8. SUMMARY AND DISCUSSION 

Since the i r initial discovery (|Toft et al 
et al.| |2005l ITrujillo et aT 



2006 



et al.||2007p the extreme properties of 



2005 D addi 
Ibtt et al. 1 12007 IZirm 



heavily debated, as galaxies with such properties were 
not expected from galaxy evolu tion models, and they do 
not exist in the local universe (jTruiillo et alp)09j |Tay~ 



lor et al. 2010 Shih & Stockton '2TJTT). Initial worries" 



that their large inferred mass densities, could be caused 
by observational biases and/or systematic uncertainties 
in their derived sizes have since been ruled out through 



higher resolution imaging (e.g. 


Cimatti et al.||2008| van 


Dokkum et al.||2008| |Toft et a 


. 2009; Damjanov et al. 


2011 


Szomoru et al. 


|201 ip . With new, more sensitive 


JN1K spectrographs li 


ke X-shooter, we are now entering 



an exciting era where it is becoming possible to detect 
absorption lines at sufficient S/N and resolution to con- 
strain their relative strengths, which are sensitive to the 
age and metallicity of the stellar populations, but not to 
dust content, and their broadening, from which we can 
derive the velocity dispersion and dynamical mass. 

In this paper we presented the first spectrum of a mas- 
sive evolved z ~ 2 galaxy of sufficient S/N and reso- 
lution to perform detailed absorption line spectroscopy. 
From the absorption lines we derive a velocity disper- 
sion of (j * = 318 ± 53km s _1 and a dynamical mass 
of log(M(j 3/ „/M0)=11.51±O.ll. Simultaneous fits to ab- 
sorption line indices and the full restframe UV-optical 
spectrum allowed us, for the first time at z <~ 2, to break 
the degeneracies between age, dust and metallicity, re- 
sulting in more realistic error bars on all derived quanti- 
ties, including the stellar mass. 



+0.2 
-0.41i 



• We derived a metallicty of \og{Z tt /Z & ) = 0.02 
which is comparable to that of similar mass local 
galaxies. 

• We accurately determined the mass/luminosity 
weighted stellar age (to within 0.1 dex), which im- 
plies that the majority of the stars formed in a 
burst at Zf > 3.3^Q'g which continued for at least 
another 0.4 Gyr (to z ~ 2.7). The derived mean 
stellar age is naturally much younger than in local 
galaxies, but if the galaxy evolves passively to the 
present day, its stellar age will be similar to local 
galaxies of similar mass. 

• We independently confirmed the galaxies quies- 
cent nature by deriving robust upper limits on the 
specific starformation rate from the slope of the 
UV continuum and upper limits on emission lines: 
sSFR < 4.4 x 10 -11 yr -1 . The shape of the con- 
tinuum and the presence of strong Balmer absorp- 
tion lines resembles a mix of A and F stars, typical 
of post star burst galaxies and consistent with the 
galaxy being quiescent, but with a significant star- 
burst relatively recently. 

We compiled 3 additional spectroscopically confirmed 
massive quiescent z ~ 2 galaxies with measured veloc- 
ity dispersion and dynamical masses from the literature, 
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allowing for a homogeneous study of scaling relations 
between structural and dynamical properties of evolved 
galaxies at this high redshift. 

Two of the galaxies are also post starburst galaxies 
with similar ages (< lGyr), brightnesses (Lx=15-23L},), 
velocity dispersions, and sizes to the galaxy analysed in 
this paper, while one is less lumnious (Lj<-=7LJj-) and 
much more compact. This galaxy has a spectrum char- 
acteristic of more evolved stellar populations, with less 
prominent Balmer lines, but stronger Ca H an K absorp- 
tion lines, and an estimated age of > 1.5 Gyr. The z ~ 2 
post starburst galaxies are not as compact as the typical 
quiescent z ~ 2 galaxies studied previously in photo- 
metric samples. This may be a selection effect, as they 
were selected for spectroscopy due to their brightness and 
therefore likely to be biased towards younger, less evolved 
z ~ 2 quiescent galaxies, where the brightness have been 
boosted by a starburst in the not so distant past. The 
evolved SEED MUSYC 1252-0 is an outlier in this pa- 
per, but may be more representative of the quiescent 
z ~ 2 galaxy population, with a fainter K-band mag- 
nitude, older stellar population and higher stellar mass 
density. The sample allowed for: 

• Studying the dynamical mass-size relation of qui- 
escent galaxies at z ~ 2, which is independent 
on systematic uncertainties and degeneracies intro- 
duced by photometric redshifts and stellar popula- 
tion synthesis modeling of broad band photometry. 

• The offset of the dynamical mass-size relation of the 
z = 2 post starburst galaxies from the local early 
type relation is smaller (a factor of 2.5) than the 
offset of the stellar mass-size relation (a factor of 
3.3), and its scatter appears to be smaller (though 
this is hard to quantify, given the small sample). 

• The smaller offsets are caused by a difference in 
stellar to dynamical mass ratios of the SEEDs and 
the local early type comparison sample. Interpret- 
ing this as a difference in the central dark matter 
fraction, the latter on average have a dark matter 
fraction of 0.46 ± 0.23 while the SEEDs on average 
have a dark matter fraction of 0.18±0.20 consistent 
with being completely baryon dominated. 



We also compared the z 
burst galaxies: 



2 galaxies to local post star- 



• The local post starburst stellar mass-size relation is 
similar to the local early type relation, but shifted 
along the relation to higher masses and larger sizes. 
The local post starburst population show a shift be- 
tween the stellar and dynamical mass-size relation, 
similar to the z ~ 2 post starburst galaxies, indi- 
cating a similarly low central dark matter fraction. 

We compared the observed structural and dynamical 
properties of our z ~ 2 sam ple to predictions from the 
cosmological simulations of Oser et al. ( 2012[ ) for the 
size, velocity dispersion and dark matter fraction evo- 
lution of massive galaxies between z = 2 and z = Q. 
The simulations do a qualitatively good job reproducing 
the difference in the stellar mass-size relations, and the 
size-velocity dispersion relations, and are also able to re- 
produce some (around 50%) of the observed difference in 
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Fig. 7. — Dynamical mass density versus the mean stellar for- 
mation redshift for the z ~ 2 galaxies. The formation redshift is 
calculated from the stellar ages, derived from the stellar population 
synthesis fits to the spectra. The error bars on the mean formation 
redshifts of the three galaxies from the literature are tentative, as 
they are not accurately derived in the papers (in one case miss- 
ing). This small sample is consistent with a picture where galaxies 
formed at higher redshift are denser. 



dark matter fraction between z=2 and 0, but not fully, 
resulting in an overprediction of the observed evolution 
in the dynamical mass-size relation. To further compare 
the structural properties of the z ~ 2 sample to local 
galaxies, we constructed: 

• The first estimate of the fundamental plane at this 
high redshift, which is offset from the local fun- 
damental plane by a larger amount than can be 
explained by passive evolution of the stellar pop- 
ulations, in other words they have to go through 
structural evolution to evolve onto the fundamen- 
tal plane observed at lower redshift. 

The structural evolution implied by the cosmological 
simulations brings the z — 2 galaxies closer to the local 
FP, but still significantly offset. This suggests that ad- 
ditional, possibly non-homologous structural evolution is 
needed. A possible mechanism responsible for this evo- 
lution is major mergers: Observations show that most 
massive galaxies will un dergo at least on e major merger 
between z=3 and (e.g. |Man et al.|2012 ), and while this 
may not be an efficie nt mechanism tor size growth (e.g 
Bezanson et al.|[2009 ), it could increase the central dark 



mat ter frac tion (th rough redist ribution of the dark mat- 
ter, Boylan-Kolchin et al.|2005 ) and provide some of the 
structural evolution needed (e.g. changing the structure 
from disk like to bulge like). 

It has been suggested that SEEDs are more compact 
than galaxies of similar mass assembled later in the his- 
tory of the Universe, because the Universe was much 
denser and m ore gas rich at earlier times 



e.g. 



Khochfar 



|fc Silkp"006| . If this is the case we may expect the mass 
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density of SEEDs to be larger f or older, more ev olved 
examples than for younger ones. Zirm et aT] (2011 1 sug- 
gested that the peak stellar density may scale with their 
formation redshift while the ratio of high-density to low- 
density components would be related to the number of 
(dry) mergers they have undergone. In Figure [8] we plot 
the central (dynamical) mass density of the z ~ 2 galax- 
ies versus their estimated formation redshift. Based on 
these four galaxies alone it is not possible to determine 
if there is a relation, but we note that the post starburst 
SEEDs have similar mass density and formation redshift, 
which are both significantly smaller than in the evolved 
SEED, consistent with this scenario. The main uncer- 
tainty in the plot is the lack of accurate ages/formation 
redshifts with realistic error bars, for the three galaxies 
from the literature. The papers describing the observa- 
tions only give rough ranges of possible ages for two of 
them, while for one (NMBS-C7447) no errors on the age 
is quoted. In a recent wor k on a photometric sa mple of 



quiescent galaxies at z ~ 2, Szomoru et al. (2011 ) did not 



find a relation between compactness and rest frame U-V 
color, which was used as a proxy for stellar age. However 
we can not rule out that this may be due to inaccuracy of 
the adopted age proxy. As shown in this paper, careful 
analysis of the absorption lines strength and the contin- 
uum emission in X-shooter spectra allows for much more 
accurate age estimates, making it possible to observa- 
tionally test for such a relation. SEED galaxies are ex- 
cellent targets for this test. They are quiescent, with star 
formation history well fitted by models with relatively 
well defined formation redshifts and they have regular 



symmetric morphologies, making it relatively simple to 
accurately determine their mass density. 

In this paper we have presented homogeneous con- 
straints on the z ~ 2 dynamical mass-size relation, the 
first constraints on the fundamental plane and a possi- 
ble relation between stellar age and compactness. Rela- 
tions that have until now only been possible study out 
to z ~ 1.3. In the immediate future it will be possible to 
put such and related studies on secure statistical footing 
as the samples of z ~ 2 — 3 quiescent galaxies with ab- 
sorption lines spectroscopy increases. This will be a large 
leap forward in our understanding of galaxy formation, 
as this is the era where these and other scaling relations 
we observe at lower redshift are believed to form. 
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APPENDIX 

STELLAR POPULATIONS: COMPARISON BETWEEN DIFFERENT OBSERVATIONAL CONSTRAINTS 

As discussed in Section[4j we derive physical parameter estimates following a Bayesian approach in which we compare 
the observed stellar absorption features, the rest-frame UV and optical color with the predictions of a large library 
of stochastic SFHs. Taking advantage of the large wavelength coverage of the X-Shooter spectrum, in this Appendix 
we explore the sensitivity of the derived parameters and associated uncertainties to the observational constraints 
adopted. In particular we repeat the analysis by using either only one color (rest-frame optical or UV) in addition to 
the absorption indices (entries 2 and 3 respectively in Table [6]) or the pixel-by-pixel flux over the rest-frame UV and 
optical range (up to 5900 A; entry 4 in Table [6]). Finally, we compare with the results obtained by fitting only the 
absorption indices to dust-free models and estimating the dust attenuation from the difference between the observed 
J-H color and the color of the redshifted, dust-free model spectra (thus effectively assuming a single-screen dust 
distribution and an attenuation law A\ oc A -0 7 , and excluding models that would predict a negative dust attenuation; 
entry 5 in Table [ 6]). T he la tter approach i s simi lar to the one adopted in the analysis of SDSS optical spectra by 
Kauffmann et al.p[2003| ) and |Gallazzi et al.| ( |2005[ ). 

In Fig. [8] we show the full one- and two- dimensional Probability Density Functions of the physical parameters of 
interests (T7e. marginalized over all the parameters except the ones shown in each panel) as obtained with the different 
observational constraints. The results are summarized in terms of median and percentiles of the one-dimensional PDFs 
in Table |H 

We note the following: 

- The luminosity- and the mass- weighted age estimates agree well within the typical uncertainty of ~0.1 dex, although 
the pixel-by-pixel fit tends to provide slightly older ages than the fits to the absorption indices alone or in combination 
with broad-band information; 

- There is a mild age-dust degeneracy in the sense that older ages are associated with lower dust attenuation. Instead, 
there is no clear age-metallicity degeneracy, partly because of the different sensitivity of the absortion indices to the 
two parameters, and pa rtly because of the la rger uncertainties in metallicity with respect to age as expected for young 
stellar populations (e.g. |Gallazzi et aL]|2005[ ). In general degeneracies between parameters are attenuated when fitting 
individual absortion features rather than the full SED; 

- The different slope in the A(FUV)-A(V) PDFs is a consequence of the different dust geometry assumed; 

- All the fits indicate a rather large attenuation by dust, although in general fits to the rest- frame optical wavelength 
range tend to predict lower overall dust attenuation, and in particular under-estimate the attenuation in the UV. The 
need to include constraints from the UV in order to correctly predict the dust attenuation is further illustrated in 
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Fig. 8. — One- and two-dimensional marginalized Probability Density Functions for the physical parameters estimated using different 
observational constraints. The filled red histograms and contours show the results of our reference fit to the absorption indices, UV and 
optical colors. These are compared to the constraints obtained by fitting the absorption indices only (blue histograms and contours; fit 5 
in Table[6]l or the pixel-by-pixel spectrum (green histograms and contours, fit 4 in Table[6]l. The contours enclose the 68%, 85% and 99% 
confidence levels. For clarity we omit the results obtained by fitting either the UV or the optical color in addition to the absorption indices, 
noting that they are bracketed by the PDFs shown here. 

Fig. [9] where the bestfit models to the five different sets of observational constraints are over-plotted on the observed 
spectrum: the best fits to the optical over-predict the flux in the UV. 

Overall there is good agreement among the various fits, in particular on stellar age. We note that our default fit 
combines the sensitivity of UV and optical colors to dust attenuation and the ability of stellar absorption features to 
alleviate parameter degeneracies. 

FITTING X-SHOOTER TEMPLATE STELLAR SPECTRA 



As an independent test of the velocity dispersion derived with pPXF we fit stellar template spectra co nvolved with 
Gaussians to the galaxy spectrum to estimate the velocity dispersion through a chi-sqared analysis (see Barth et al. 
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TABLE 6 

Stellar population parameter estimates and uncertainties obtained by fitting 
different observational constraints. the results adopted in the paper are those 

obtained by simultaneously fitting three stellar absorption features, the uv 
flux ratio and the observer-frame j-h color (bold face). the other fits explore 
the effect of changing the observational constraints (see text). 



Method 


log(t r /yr) 
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log(Z»/Z Q ) 


log(M*/M ) 


(1) indices, UV and J-H color 




9.08tg;iS 


0.0212.41 


11 S7+ 013 
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(2) indices, J-H 

(3) indices, UV flux ratio 

(4) full SED 

(5) indices, color excess 
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TV 




Apuv 
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(1) indices, UV and J-H color 
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(2) indices, J-H 

(3) indices, UV flux ratio 

(4) full SED 

(5) indices, color excess 


1 oy + U.Sl 
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2 04+ 101 
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TABLE 7 

Stellar templates used in velocity dispersion 
fitting. The best fit velocity dispersion AND 68 PER 

CENT CONFIDENCE INTERVAL IS LISTED IN THE LAST TWO 

COLUMNS. All fits were DONE in the region between 
3785 AND 5000 A (restframe) as explained in the text 



Template 


Spectral type mix (%) 


£7 


68% CI 




AV 


F5V 


[km s -1 ] 


km s _1 


1 


100 





300 


210,380 


2 





100 


290 


240,350 


■1 


75 


25 


270 


200,340 


5 


50 


50 


260 


200,330 


6 


25 


75 


270 


220,330 



2002| |Wold et al.| [2007). For stellar templates we use stars of five different spectral types (AV, F5V, G81V, K0III, 
K0V), observed with X-shooter at the same slit width, for a different program (084.B-035A, PI: Hjorth). In order to 
estimate the velocity dispersion of the galaxy we convolved these stellar spectra with Gaussians and fitted them to the 
Balmer line region of the galaxy spectrum 

The rest-frame wavelength range of the NIR arm in the galaxy spectrum corresponds approximately to the wavelength 
range of the UVB arm in the stellar spectra. When we do the spectral fitting, we therefore compare the stellar UVB to 
the NIR galaxy spectrum. The stellar spectra were reduced with a spectral sampling of 0.2 A/pix in the UVB arm, and 
0.5 A/pix in the NIR arm. The spectral sampling in the NIR arm of the galaxy spectrum corresponds to 0.16 A/pix 
in the galaxy rest-frame. As this is very similar to the sampling of the UVB stellar spectra, and because of the lower 
quality of the galaxy spectrum compared to the stellar spectra, we chose to use the spectra as they are without further 
manipulation and resampling. Both the galaxy spectrum and the stellar spectra were normalized by their continua 
before fitting. For the stellar spectra we extracted the continuum by marking anchor points in the spectra and doing a 
spline fit with varying degrees of tension. Due to the lower S /N we used a different approach for the galaxy spectrum, 
and extracted the continuum shape from the best fitting model of the galaxy (see Section [4]) with sampling every 50 
A. We experimented with fitting different regions of the spectrum to maximize the S/N in the velocity dispersion 
measurement, including as many absorption lines as possible, while down-weighing regions of the spectrum with low 
S/N or strong skyline residuals. We obtained the best results by fitting the region between rest-frame wavelength 3785 
A and 5000 A where most of the strong absorption features are found. This region covers the region bluewards of 
H(3. The Mgl triplet at «5175 A redwards of H/3 and the Na doublet (5889.95, 5895.92 A) are not detected in the 
spectrum, hence we concentrated on estimating the velocity dispersion from the Balmer line region. The AV and F5V 
stars (which have prominent Balmer lines) provide good fits with well constrained minima in the x 2 distribution. The 
other stars provide significantly worse fits. The KOV has a minimum consistent with the best fitting values of the K0V 
and F5V star, but with a much flatter \ 2 distribution, while the G8IV and K0III do not produce a minimum, due to 
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(5) indices, color excess 




4000 



5000 



6000 



7000 




2000 



2200 



2400 



2600 



2800 



3000 



3200 




200 



300 



400 



1500 
X [A] 



600 



700 



800 



Fig. 9. — X-Shooter spectrum at medium-resolution (3A/bin in the NIR arm, 16A per bin in the UVB and VIS arms) over the full 
wavelength range covered by the three arms. Overplotted are the bestfit model spectra obtained by constraining 1) three stellar absorption 
indices, one rest-frame optical color and the UV flux ratio (red, our default fit), 2) three stellar absorption indices and one rest-frame 
optical color (magenta), 3) three stellar absorption indices and the UV flux ratio (cyan), 4) the full spectrum pixel-by-pixel up 5900A 
(green), and 5) three stellar absorption indices and a cut on color excess (blue). Fits obtained by constraining the optical wavelength range 
only under-predict the dust attenuation in the UV and hence over-predict the UV flux. In the NIR spectrum, we have masked out areas 
of high background contamination. 

too weak Balmer lines, and too strong G-band and CaH+K lines. We also fit different linear combinations of the AV 
and F5V stars. These fits produce consistent best fitting values, but with slightly smaller errorbars. In Figure [lO] and 
Table [7] we show the results from the best fits. From this analysis the best fitting velocity dispersion is in the range 
260 — 300 km s _1 with a typical (random) uncertainty of 50km s -1 , in agreement with the results from the pPXF fits. 
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Fig. 10. — Fits to the galaxy spectrum with the A and F star, and the best fitting mixture. Regions excluded in the fit are shaded in 
grey. The galaxy spectrum is smoothed to a resolution of 10 A for displaying purposes. 
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